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Understanding and predicting 'decoherence' is impor- 
tant in future applications involving nanonictric devices, 
fast switching times, clock synchronization, superdense 
coding, quantum computation, teleportation, quantum 
cryptography, etc. where entanglement structure and dy- 
namics play a key role . In the last three decades it has 
also been central in exploring possible limitations to the 
validity of standard unitary quantum mechanics (QM), 
by studying a variety of linear and nonlinear extensions 
that have been advocated by several authors on a vari- 
ety of conceptual grounds |2l|. It has been suggested ^ 
that long-baseline neutrino oscillation experiments may 
provide means of testing the existence of spontaneous 
decoherence at the microscopic level and the validity 
of linear and nonlinear extensions of the Schrodinger- 
von Neumann equation of motion of QM, thus prompting 
a renewed interest on such extensions |^ U, y| . 

The aim of this Letter is to list the main conditions 
that must be imposed and checked on linear and nonlin- 
ear extensions of QM which assume an augmented set of 
true quantum states described by state operators p with- 
out the restriction = p. The reasoning and framework 
proposed here should provide useful guidance also to cur- 
rent efforts to define general measures of entanglement 

H 

The conditions proposed here form a very restrictive 
set. Yet, at least one possible extension has been proved 
to satisfy them all'r'l, with mathematics that has been 
partially rediscovered recently by researchers in different 
contexts and fields 0, • 

1. Causality. Forward and backward in time. 

We consider the set V of all linear, hermitian, 
nonnegative-definite, unit-trace operators p (without the 
restriction p^ — p) on the standard QM Hilbert space 
Ti associated with a strictly isolated system[l(|. Every 
solution of the equation of motion, i.e., every trajectory 
u{t, p) passing at time t — through state p in V, should 
lie entirely in V for all times t, —oo < t < +oo. This 
strong causality condition is nontrivial and demanding 
both from the conceptual and the technical mathemati- 
cal points of view. 

2. Conservation of energy and other invariants. 



The value of the energy functional e(p) = Tr{pH), where 
H is the standard QM Hamiltonian operator associated 
with the isolated system [H ^ H{t)], must remain in- 
variant along every trajectory. If H is the Fock space 
of an isolated system consisting of M types of elemen- 
tary constituents (e.g., atoms and molecules if chemical 
and nuclear reactions are inhibited; or atomic nuclei and 
electrons for modelling chemical reactions) each with a 
number operator Ni {[H, Ni] = and [N,,Nj] = 0), then 
also the value of each number-of-constituents functional 
n.i(p) — Tr(piVi) must remain invariant along every tra- 
jectory. Depending on the type of system, there may be 
other time-invariant functionals, e.g., the total momen- 
tum components Pj{p) = Tr(pPj), with j = x,y,z, for 
a free particle (in which case Galileian invariance must 
also be verified, for [H,Pj] = and [P^,Pj] = 0). In 
what follows, we denote by gi{p) = Tv{pGi) the set of 
non-Hamiltonian time-invariant functionals, if any, with 
[H,G,] = and [Gi,Gj] = (clearly, H and the G,'s 
have a common eigenbasis that we denote by {IV'f)})- 
3. Standard QM unitary evolution oi p^ — p states. 
Unitary time evolution of the states of QM according to 
the Schrodinger equation of motion must be compatible 
with the more general dynamical law. These trajecto- 
ries, passing through any state p such that p^ = p and 
entirely contained in the state domain of QM, must be 
solutions also of the extended dynamical law. Because 
the states of QM are extreme points of the state domain 
V, the trajectories of QM must be boundary solutions 
(limit cycles) of the extended dynamical law. 

In general, any extended dynamical equation may be 
written in the form 

^^-^[H,p]+Dm (1) 
withi^M -DAf(p,i^,G„...) , (2) 

where operator D]^J represents the dissipative part of the 
equation of motion and may depend linearly and/or non- 
linearly (through superoperator Da/) on the state oper- 
ator p, on the Hamiltonian H, on the linear operators 
Gi associated with the other time invariants (if any) , as 
well as on the structure and the number M of elemen- 
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tary constituents of the system. Like the Schrodinger- 
von Neumann term, also the dissipative term should not 
be responsible for rates of change of any of the invariant 
functionals Tr(p), e(p), gi{p) and, therefore, 

TtDm = TtDmH = TrDMGi = . (3) 

If the complete dynamics preserves the feature of 
uniqueness of solutions throughout the state domain 
V, then pure states can only evolve according to 
the Schrodinger equation of motion and, therefore, 
Dm(p, H,Gi, . . .) — when = p. This feature that 
may be responsible for hiding the presence of deviations 
from QM in experiments where the isolated system is 
prepared in a pure state. It also implies that no trajec- 
tory can enter or leave the state domain of QM. Thus, 
by continuity, there must be trajectories that approach 
indefinitely these boundary solutions (of course, this can 
only happen backward in time, as t — > — oo, for other- 
wise the entropy of the isolated system would decrease in 
forward time) . 

4. Conservation of effective Hilbert space dimen- 
sionality. 

Unitary dynamics [Eq. Q with Dm — 0] would main- 
tain unchanged all the eigenvalues of p and therefore can- 
not satisfy Condition 5 below 0|. Instead, we only re- 
quire that the dynamical law maintains zero the initially 
zero eigenvalues of p and, therefore, conserves the car- 
dinality of the set of zero eigenvalues, dimKer(/9). In 
other words, if the isolated system is prepared in a state 
that does not 'occupy' the eigenvector li/ji) of H (and the 
Gi's), i.e., if p{0)\ipe) = (so that is also an eigen- 
vector of p corresponding to a zero eigenvalue) , then such 
energy eigenvector remains 'unoccupied' at all times, i.e., 

pit)\A) = 0. 

This condition preserves an important feature that al- 
lows remarkable model simplifications within QM: the 
dynamics is fully equivalent to that of a model system 
with Hilbert space 7i' (a subspace of H) defined by the 
linear span of all the \tpe)'s such that p{t)\ipe) 7^ at 
some time t (and, hence, by our condition, at all times). 
The relevant operators X' on H' {p' , H' , ...) are 
defined from the original X onH {p, H, Gi, . . .) so that 
{ak\X'\ae) — {ak\X\ae) with \ak) any basis oiH'. 

It is also consistent with recent experimental tests 
that rule out, for pure states, deviations from linear and 
unitary dynamics and confirm that initially unoccupied 
eigenstates cannot spontaneously become occupied. This 
fact adds nontrivial experimental and conceptual difh- 
culty to the problem of designing a fundamental test of 
QM, capable of ascertaining whether decoherence orig- 
inates from uncontrolled interactions with the environ- 
ment due to the practical impossibility of obtaining strict 
isolation, or else it is a more fundamental intrinsic feature 
of microscopic dynamics requiring an extension of QM. 
In the latter case, this condition will preserve within the 



extended theory the exact validity of all the remarkable 
successes of QM. 

5. Entropy nondecrease. Irreversibility. The prin- 
ciple of nondecrease of entropy |0] for an isolated system 
must be satisfied, i.e., the rate of change of the entropy 
functional — A:BTr(plnp) must be nonnegative along ev- 
ery trajectory, —kjf£T:[u{t,p)hiu{t,p)] > — A:BTr(/3lnp). 

6. Stability and uniqueness of the thermodynamic 
equilibrium states. Second law. A state operator p 
of the isolated system represents an equilibrium state if 
dp/dt = 0. For each given set (e, g) of feasible values of 
the energy functional e{p) and the other time- invariant 
functionals gi{p), if any, among all the equilibrium states 
that the dynamical law may admit there must be one and 
only one which is globally sto&fe|l3|. 

This stable equilibrium state must be that of equilib- 
rium thermodynamics and, therefore, of the form 

^ exp[-/3(e, g)^ + X], i^zje, g)G»] , . 

Trexph/3(e~,g)i/ + E.'^.(e,g)G.] ' ^' 

where Gi are defined above. Of course, states given by 
Eq. Q are solutions of the constrained maximization 
problem 

max — fcBTr(plnp) subject to (5a) 
Trip) = 1, TiipH) = g, Tr(pG,) = g,, and p > (Bb) 

and reduce to the canonical equilibrium states pc ~ 
exp(-/3i?)/Trexp(-/3i?) when Gi = gj for all Vs (with 
gi scalars and / the identity on 7i), and to the mi- 
crocanonical state pc — 1/ dimH if also H ^ el (and 
dimTi < 00). 

As discussed in Ref. 0| , the entropy functional is not a 
Lyapunov function, even if, in a strict sense that depends 
on the continuity and the conditional stability of states 
Pc it does provide a criterion for the local stability of 
these states. In addition to this, the second law requires 
however that no other equilibrium state of the dynamical 
law be be globally stable[l3llli|. 

Consider the noteworthy family of states 

_ BeM-P{e.m + l:^^^{^.k)G^]B 
^"'^ TrSexp[-/3(e~,g)i? + E.i^.(e,g)G.] ' ^ ' 

where B is any given idempotent operator B^ = B. 
This family, which includes pure states [Tr(i3) = 1], 
maximizes the entropy [Prob. subject to the addi- 
tional constraint p = BpB for the given B. All eigen- 
values of pnd must remain invariant (otherwise the en- 
tropy would decrease) and the state is equilibrium if 
[B,H] — or otherwise it evolves unitarily (limit cycle) 
with B{t) = exp{-iHt/h)B{0) exp{iHt/H). They have a 
thermal-like distribution (positive and negative tempera- 
tures) over a finite number [Tr(i3)] of 'occupied' eigenvec- 
tors. Because entropy cannot decrease and — fcBTr(plnp) 
is an S'-function, they are conditionally locally stable 
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equilibrium states or limit cycles [l5j- For them not to 
be globally stable, as required by the second law, it suf- 
fices that the extended dynamics imply that at least one 
state of equal energy and other invariants (not necessar- 
ily neighboring nor with the same kernel) evolves towards 
higher entropy than /9„d- 

7. Non-interacting subsystems. Separate energy 
conservation. For an isolated system composed of 
two distinguishable subsystems A and B with associated 
Hilbert spaces Ti,^ and , so that the Hilbert space 
of the system is H = H^'S)'H^ , if the two subsystems 
are non-interacting, i.e., the Hamiltonian operator H — 
Ha®Ib +Ia'^Hb, then the functionals Tr[{HA®lB)p] = 
Tia{HaPa) and Tr[(/A®i?s)p] = ^j:b{HbPb) represent 
the energies of the two subsystems and must remain in- 
variant along every trajectory, even if the states of A 
and B are correlated, i.e., even if p 7^ Pa®Pb- Of course, 
PA = ^t:b{p), Pb = Tr^(/9), Tr^ denotes the partial trace 
over Ti^ and Tr^ the partial trace over Ti^. 

8. Independent states. Weak separability. Sepa- 
rate entropy nondecrease. Two distinguishable sub- 
systems A and B are in independent states if the state 
operator p — pa®Pb- For any given p, let us define 
the idempotent operator B obtained from p by substi- 
tuting in its spectral expansion each nonzero eigenvalue 
with unity [I4I and the entropy operator S = —k^Blnp 
(always well-defined). For independent states, S = 
Sa^Ib + Ia(®Sb = -k^[BA\ripA®lB + Ia^BbI^pb]- 
For permanently non-interacting subsystems, every tra- 
jectory passing through a state in which the subsystems 
are in independent states must proceed through indepen- 
dent states along the entire trajectory, i.e., when two un- 
correlated systems do not interact with each other, each 
must evolve in time independently of the other. 

In addition, if at some instant of time two subsystems 
A and B, not necessarily non-interacting, are in inde- 
pendent states, then the instantaneous rates of change of 
the subsystem's entropy functionals —ksTT:{pA In pa) and 
—k^Tv[pB^T^ Pb) must both be nondecreasing in time. 

9. Correlations, entanglement and locality. 
Strong separability. Two non-interacting subsystems 
A and B initially in correlated and/or entangled states 
(possibly due to a previous interaction that has then been 
turned off) should in general proceed in time towards less 
correlated and entangled states. In any case, in order 
for the dynamics not to generate locality problems, i.e., 
faster-than-light communication between noninteracting 
subsystems (even if in entangled or correlated states), 
entanglement and correlations must not increase in the 
absence of interactions. In other words, when subsystem 
A is not interacting with subsystem B, it should never be 
possible to infiuence the local observables of A by acting 
only on the interactions within B, such as switching on 
and off parameters or measurement devices within B. 

This however does not mean that existing entangle- 



ment and/or correlations between A and B established 
by past interactions should have no infiuence whatso- 
ever on the time evolution of the local observables of 
either A or B. In particular, there is no physical rea- 
son to request that two different states p and p' such 
that p'j^ — PA should evolve with identical local dynamics 
(d/9^/dt — d/9A/dt) whenever A does not interact with _B, 
even if entanglement and/or correlations in state p differ 
from those in state p' . Rather, the two local evolutions 
should be different until spontaneous decoherence (if any) 
will have fully erased memory of the entanglement and 
the correlations established by the past interactions now 
turned off. In fact, this may be a possible experimental 
scheme to detect spontaneous decoherence. 

Compatibility with the predictions of QM about the 
generation of quantum entanglement between interact- 
ing subsystems that emerge through the Schrodinger- 
von Neumann term —i[H,p]/h of Eq. requires that 
the dissipative term Dm may entail (spontaneous) loss 
of entanglement and loss of correlations between subsys- 
tems, but should not be able to create them. 
10. Onsager reciprocity. First, we introduce a par- 
ticularly useful representation of general nonequilibrium 
states 7]. Given any state p on H, we define the effective 
Hilbert space H' as above, and choose a set of operators 
{/', X[,X2, ■ ■ ■} spanning the linear space ChiH') of lin- 
ear hermitian operators on 7i'; the corresponding state 
p' on Ti' has no zero eigenvalues, so that S — —kj^Blnp 
becomes S" = — feglnp' on W, which can be written as 
S' = fof + ZjjX'j because it belongs to ChiH'). Thus, 



exp(-S,-/j^i/fcB) 
Trexp(-E,/,Xj/fcB 



(7) 



where /o = fee lnTrcxp(— /^Xj/fce). Similarly, we 
can also write S'^ — —kg In p'^ — foel' + T,jfjeX'j, for the 
target maximum-entropy equilibrium state on 



P'eip') 



Trexp(-/3i?' + E,.i.fcG'fe) ' 



(8) 



where (3 and are such that e(p'^) — e(p') and 
M) = 9^{p'), so that Tr(p'lnp^) = Tr(p^lnp^) and 
Tr[(dp'/dt)iS'g] =0. As a result, the following relations 
can be easily proved, 



S{p') - S{p',{p')) =fo- foe + E,(/. - /.e)x.(p') , (9) 

(10) 



d[s{p) - S{p'^{p))] _ f _ r 

dx,{p') - ' 



-l^iJ^ Dt - l^iUt lie) 

(A^'A^) = E., f J, (AX^AXj) , (12) 

where Dxi{p')/Dt = Tr{DMX^) denotes the dissipa- 
tive rate of change of the linear mean-value functional 
x,{p') = Trip'Xl), {AS' AS') = Tr[p'(-fc3 Inp')'] - 
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s{pr, (AX;AXj) = ^Tr[{p{Xl,X'^}] - x,{p')x,{p'). 
When the system is in state p' , we interpret (AX-AXj) 
as the codispersion (covariance) of simultaneous measure- 
ments of observables X^ and Xj, (AX^AX^) as the dis- 
persion (or fluctuations) of observable X^ and {AS' AS') 
the entropy fluctuations. 

By Eq. (|10|l . we may also interpret fi — fie as the gener- 
alized affinity or force conjugated with the mean value of 
the linear observable Xi. In order to recover Onsager's 
theory, we may impose that (at least in the vicinity of 
state Pg) the extended dynamics be such that the dissipa- 
tive rates be linearly related to the generalized affinities 
through generalized-conductivity functionals, i.e., 

= J^i., (p', i?', G'„X',, ...)(/, - , (13) 

3 

where the i^j's may be nonlinear functionals of p' (pos- 
sibly to be approximated with their values at p^, in its 
vicinity), but should form a symmetric (H — > — H, if 
H' depends on an external magnetic field) non-negative 
definite matrix, so that the rate of entropy production 
results in a quadratic form J2zjifi ~ fie)Lij{fj - fje)- 
Moreover, the iy's should be linearly interrelated with 
the matrix of codispersions {AX^AX'^), in order to re- 
cover also Callen's fluctuation-dissipation theorem. 
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